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Abstract 



To understand the short versus long-distance final-state interactions, we 
have performed a detailed amplitude analysis for the two-body decay J /ip — > 
1 _ 0~. The current data favor a large relative phase nearly 90° between the 
three- gluon and the one-photon decay amplitudes. The source of the phase 
is apparently the long-distance final-state interaction. Nothing anomalous is 
found in the magnitudes of the three-gluon and one-photon decay amplitudes. 
We discuss the implications of this large phase in the weak decay of heavy 
particles. 
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I. INTRODUCTION 



Though the final-state interaction phases are important to observability of CP violating 
decays, it is very difficult to compute them or to extract them from data. Only the short- 
distance contribution has been computed in the quark-gluon picture Any attempt 
to estimate the long-distance contribution in the hadron picture has so far been limited to 
elastic or quasi-elastic rescattering, which is presumably only a small portion of the long- 
distance effect, particularly in the heavy particle decays like B-decay. Some argue that 
the long-distance contribution is negligible at least in some decay modes ||, while others 
identify a specific long-distance contribution and show that it is actually much larger than 
the short- distance effect . Some entertain the idea that many long-distance contributions 
might average out to a small effect after being summed up. 

In this paper we try to test whether the short- distance final-state interaction dominates 
over the long-distance one in the J/i/j decay or not. Though the J/ip decay proceeds with 
strong and electromagnetic interactions, its narrow width allows us to treat the decay just 
like weak decays, namely, a short- distance decay followed by long-distance rescattering. In 
the quark-gluon picture J/ if) decays either directly into three gluons or into a quark and an 
antiquark through one photon. Both processes acquire a short-distance QCD rescattering 
phase of 0{o, s /'k). If this is the dominant source of the final-state interaction phase, the 
relative phase ought to be very small for the amplitudes of all decay modes. On the other 
hand, if long-distance processes are important to generating the phase, the decay amplitudes 
can have large relative phases to each other. Since there are sufficient data on the decay 
J/ip — > 1~0~, we are able to perform an amplitude analysis and extract the relative phase 
between the three-gluon decay and the one-photon decay amplitudes. 

In Section II we attempt a detailed numerical analysis of the decay amplitudes in the 
framework of the broken flavor SU(3) symmetry. We include all first-order symmetry break- 
ings and some of second order effects. The result of our analysis shows that the three-gluon 
and one-photon amplitudes have a large relative phase of rescattering, nearly 90° off phase 
to each other. It indicates that the major source of the final-state interaction phases is in 
the long-distance hadronic rescattering. 

The decay J/ip — > 1~0~ has been a subject of discussion in connection with the abnor- 
mally small yield of the decay ip' — ► 1~0~. A few exotic models were proposed to resolve 
this puzzle |5]||0]. They suggested that the dominant process of J/ip — > l - 0~ is not a per- 
turbative three-gluon decay. However the result of our analysis shows that the magnitude 
of the 1 = decay amplitudes is consistent with the three-gluon decay. We can also show a 
serious shortcoming of those models. We shall discuss on this point in Section III. Finally 
in Section IV, we discuss on the implications of the present analysis in the B-decay. 

II. AMPLITUDE ANALYSIS IN BROKEN FLAVOR SU(3) SYMMETRY 

Before entering our amplitude analysis, we read off one relevant information from the 
current data. The lepton-pair decay branching fraction B\ for J/ip — > 7 — > l + l~(= e + e~ + 
has been measured with a high accuracy. In contrast, the inclusive one-photon 
annihilation into hadrons for J/ip — > 7 — > qq(= uu + dd + ss) was obtained only indirectly 
from the nonresonant background cross section a(e + e~ — *> qq) interpolated to the J/ip mass. 
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This two decades old value || quoted in the Review of Particle Physics (RPP) || has a large 
uncertainty. Actually we can obtain a more accurate value for B 1 = B{J/ip — > 7 — > qq) 
by making the perturbative correction of (1 + a s /ir) to the leptonic branching fraction Bi. 
With this help of theory, we can compute the inclusive three-gluon decay branching fraction 
B g gg with B ggg = 1 — Bi — £> 7 . We thus obtain for the inclusive three-gluon and one-photon 
annihilations into hadrons 

= 5.8 ±0.3 (4.2 ±0.6), (1) 

where the number in the bracket is the experimental value quoted in RPP. Except for the 
number in Eq.(|l]), we take the numbers tabulated in RPP as the current data throughout 
this paper. 

If branching fractions of the J/ip —> 1~0~ modes scale with those of the inclusive hadron 
decay processes, we expect that the amplitude analysis should give us the ratio of the three- 
gluon to one-photon decay amplitudes into the exclusive channels 1~0~ somewhere between 
2 and 3 (~ \^x8). Later we shall often compare the results of our amplitude analysis with 
this number. 



A. Fit to J/ip -» V 9 P 8 

We first study the two-body decay modes of J/ip — > 1~0~ in which the pseudoscalar 
meson belongs to an SU(3) octet. The singlet rf will be included later with the T] — rf 
mixing. For the vector mesons, we study the singlet and the octet together as a nonet. We 
parametrize the decay amplitudes as follows: 

(1) The 1~ mesons form the ideally mixed nonet ||10|| , namely <fi = —ss. Therefore the 
SU(3) symmetric coupling is given by 

L int = atr(V 9 P 8 ), (2) 

where V 9 and P$ are represented in the 3x3 matrices. 

(2) The strong SU(3) breaking of As is included in the three-gluon decay: 

L mt = etr({V 9 ,P 8 } + T 3 3 ), (3) 

where we use T 3 3 instead of As to simplify the numerical coefficients of parametrization. The 
symmetrization of V 9 and P$ is required by charge conjugation invariance. 

(3) The one-photon annihilation amplitudes transform like A^ = (A3 + As/v / 3)/2. There- 
fore, they are parametrized as: 



J int 



a y tr({V 9 ,P 8 } + X E ). (4) 



(4) The phases of the amplitudes a, e, and a 7 are group theoretically independent. There- 
fore, we introduce two relative phases, <5 7 and S e , defined as: 

arg(a 7 a*) = <5 7 arg(ea*) = S e . (5) 
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(5) The p — uj [|TTJ mixing can be potentially important to the processes of AJ = 1. For 
instance, the process J '/if) — > lot] — > pr] interferes with the direct process J / 'ip — > 7 — > pi] and 
is counted as part of the branching to J /if) — > prj. The decay J '/if) — ► lot] can proceed through 
three gluons while the direct process is the one-photon process. Therefore, the p — u mixing 
may not be negligible in this mode. Since the 10 width is much narrower than the p width, 
the major contribution of the p — uj transition to J /if) — > pi] occurs at the to resonance peak 
of 7r + 7r~. In comparison, in the process like J/if> — > ujtt , the effect of the p — uj transition 
is less significant because the large p width suppresses it kinematically. In the processes of 
AI = 0, the p — uj transition is negligible. 

We have included the p — uj mixing with the effective transition coupling, 

Unt = mlfp^p^, (6) 

with fpu = —6.8 x 10~ 3 , as extracted from the process e + e~ — > tt + tt~ around the uj mass. 

(6) The p-wave phase space correction is made with p 3 for all decay branching fractions. 
If the flavor symmetry applies best to the dimensionless decay couplings, the phase space 
factor p 3 should be divided by some quantity having the dimension of squared mass that 
may be subject to symmetry breaking. However, this uncertainty has been incorporated 
through the As breaking of (2) above, at least, to the first order. 

In Table I, we have tabulated the parametrization of seven VgPg decay amplitudes with 
the SU(3) amplitudes, 

a, e, a 8 . (7) 

If the 1~0~ decay branchings scale more or less with the inclusive ones (cf Eq. ([[])), we 
expect |a 7 /a| ~ 0.67 in the normalization used in Table I. For the strong symmetry breaking, 
\e/a\ < 0.3 is a reasonable range. In addition, we have two relative phases as free parameters, 
S t and 5 7 . 

The best fit with these five parameters is obtained for 

a = 1, e = -0.22, a 7 = 0.34, 5 e = -22.5°, 5 7 = 80.3°, (8) 

where magnitude of the amplitude a has been normalized to unity. The fitted values to the 
data are tabulated in the third column of Table next to the observed values. The x 2 is 4.8 
for this fit. Though the ratio of a 7 /a in Eq.(|8|) is a half of the scaled value ~ 0.67, it can 
hardly be called an enhancement of the three-gluon decay. The magnitudes of a, a 7 , and 
e are in line with the expectation from the inclusive branching ratios: There is no sign of 
significant enhancement of the three-gluon processes relative to the one-photon processes. 

If we fit the data without the phases, the best \ 2 is 57. When there is no rescattering 
phase, the p — uj mixing is unimportant because the p — uj transition amplitudes are 90° 
off phase to the main amplitudes. The fitted values are listed in the fourth column (No 
phase I) of Table II. In order to make a quantitative comparison of the fits with and without 
the phases, it is more appropriate to fit the data with the same number of free parameters, 
namely five real amplitudes. We may add the amplitudes of the As breaking to a 7 as the 
second-order small quantities. We may also include breakdown of the ideal nonet coupling 
ansatz. Actually, there is a subtlety between breakdown of the nonet and the 0(A§e) cor- 
rection. In the nonet scheme, which is realized in the nonrelativistic quark model, V\ does 



4 



not form an SU(3) coupling without accompanying Vs. For instance, the term like trVg oc V\ 
is not allowed in an SU(3) coupling. When a A 8 breaking is taken into account, however, 
one may include tr(A 8 Vg) or tr(T 3 3 K,) among others. Normally it would not matter which 
of these is included since their difference has the same SU(3) property as the term of the 
symmetry limit. In the case of the nonet, the difference ~ trVg is a term forbidden by the 
nonet coupling ansatz. Therefore we must choose between them either from the observed 
Vg parameters or from some theoretical reasoning. Since the A 8 breaking is caused by the 
s-quark mass term and the orgin of the nonet is in the nonrelativistic quark model, we feel 
that T| is more appropriate than A 8 in the case of the nonet. We include the strong SU(3) 
breaking to a 7 along this line to see an outcome of the fit. 

After charge conjugation invariance is taken into account, there are three independent 
amplitudes of 0(T 3 3 e) that have different group structures from the amplitudes in Eq.(^). 
One of them (oc tr(T 3 3 Vg)) contributes only to the 07r° and (prj modes. Meanwhile a severe 
upper bound has been set on B(J/ip — > 07r°) by experiment. With this upper bound, we 
find that this SU(3) amplitude contributes no more than 6% to 4>r), smaller than normally 
expected for strong SU(3) breakings and below the level of our concern. We therefore drop 
this small amplitude and retain only the remaining two amplitudes as the 0(T 3 3 e) corrections: 

L int = e 7l [tr(V 9 A i; P 8 r 3 ) + tr(P s \ E V 9 T*)} + e j2 tr(P 8 T*)tr(V 9 \ E ). (9) 

Now we have five real amplitudes. The best fit attains \ 2 = 21 in this case. The fitted 
values are: 

a = 1, e = -0.14, a 7 = 0.30, e 7 i = -0.12, e 72 = -0.11. (10) 

The fitted branching fractions are listed in the fifth column (No phase II) of Table II. 

A simple qualitative explanation can be given as to why the best fit needs the phases. 
Refer to the Table I and the observed branching fractions listed in Table II for the following 
discussion: 

1. First of all, the significant difference between the pn and K + *K~ branchings must 
be explained by the SU(3) breaking e amplitude. The e amplitude of a right magnitude 
(|e/a| = 0.22) produces this difference. 

2. Next, we need a sizable a 7 amplitude in order to account for the utt° mode. 

3. The a 7 amplitude contributes to splitting the branching fractions of K + *K~ and 
K°*K too. However, if a 7 and a substantially interfered, this splitting would be much too 
large. To keep the splitting between K + *K~ and K *!^ small, a 7 and a must be largely off 
phase to each other. 

The relative phase between e and a turns out to be small. One interpretation for the 
smallness of this relative phase is that the main source of e is the kinematical SU(3) breaking 
due to mass splitting in the phase space and decay coupling. 

B. Including 7/ 

Once rj is included, the number of independent parameters suddenly increases since 
unlike the vector meson couplings, the singlet 0~ couplings are not related to the octet 0~ 
couplings. Complication grows further when we include the rj — rj' mixing. 
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Here we present a relatively simple sample of analysis instead of the most general one in 
order to show that need of the relative phase between the three-gluon and the one-photon 
amplitudes persists. We add all SU(3) independent amplitudes involving rf that correspond 
to a, e, and 7 of P$: 

L mt = a'tr^Px) + e'tiiV^Tl) + a^r^P^). (11) 

In addition to the p — to mixing, we include the 77 — 77' mixing, 

r] = T] 8 cos 6 P - r]i sin 9 p , r( = r] 8 sin 9 p + r] 1 cos 8 P , (12) 

where 9 P is —10° ~ —20° [|T2 1 . For the relative phase, we put a common phase 8 between all 
one-photon amplitudes and all three-gluon amplitudes {e.g., 5 e = 0). 

The best fit to the data is obtained with the following values of the parameters: 



a = 1, e = -0.18, a 7 = 0.36, (13) 
a 



' - 0.44, e' = 0.051, a' = -0.40, 



"1 " 
5 = 75.2°, 

The value of x 2 is 8.2 for fitting eleven data with seven parameters. The relative magni- 
tudes of the parameters are normal, namely, in line with the expectation from the inclusive 
branchings and strong SU(3) breakings. On the other hand, if we attempt to fit the data 
without the phase 5, the x 2 jumps to 43. In this case, sum of the one-photon branchings is 
close to that of the three-gluon branchings for VgP\. The tendency of deterioration of the 
fit without a phase persists as we have seen in the case without 77'. The fitted values of the 
branching fractions are tabulated with and without the phase in the last two columns of 
Table II. The phase is unimportant in fitting to the decay modes involving rj'. 

We conclude that as long as the currently listed data are taken at their face value, the 
three-gluon and one-photon amplitudes have a large relative phase to each other. Apart 
from this unexpected result, our amplitude analysis show that the magnitudes of all decay 
amplitudes are within the range of what we expect. 



III. IS THE DECAY J/tp -> 10 ANOMALOUS ? 

We have chosen the J/tp — » 1~0~ decay modes for study of the final-state interaction 
phases since they are the most extensively measured decay modes. No similar analysis can 
be made for other modes at present. 

Meanwhile, there was one disturbing twist related to these decay modes. That is, the 
ip' — > 1~0~ decay modes are severely suppressed in comparison with the corresponding J/ip 



modes [13]. For pn, the upper bound on the branching fractions normalized to the e + e 



branching fraction obey the inequality, 



B »''^^ <1.2xl0-»x B M^k- (") 



B(ip'^e+e-) ' B(J/^^e+e~ 

This vast difference between J/ip and ip' has stimulated many speculations on the pure QCD 
decay of J/tp and ip'. The argument goes as follows: Normally J/ip(ip') — > 1~0~ would be 
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highly suppressed by chirality mismatch of perturbative QCD. However, this suppression is 
compensated by an enhancement in the 1 = channels of J /if). The enhancement brings 
T(J/if) — > 1~0~) back to the value predicted by the perturbative three-gluon annihilation 
0. The cause of enhancement may be either a vector gluonium state nearly degenerate 
with J/if) or an hidden charm pair in the light 1~ mesons 0. T(if)' — > 1~0~) is small 
because it suffers from chirality mismatch but receives no enhancement. 

However, our amplitude analysis raises a doubt about such an explanation. We have seen 
that both the 1 = and the 1 = 1 amplitudes are as normal as we expect from the inclusive 
three-gluon and one-photon annihilations. If the observed magnitude of the 1 = amplitudes 
were actually the result of the compensation between a chirality suppression and a dynamical 
1 = enhancement, we would expect that the one-photon annihilation amplitude for u)n(I = 
1) should be suppressed by chirality without a compensating enhancement. If their models 
are correct, we can read off the chirality suppression factor from Eq . flTil) . With the chirality 
suppression, B(J/ip — > un ) would have to be 

B{J/if) -> uir°) « (chirality suppression) x ^^Z^l 1^1 . B{J/if) -> pvr) (15) 
1 

in those models. The data violate this inequality by two orders of magnitude. Therefore, the 
origin of the relative suppression of if)' — > 1~0~ to J/if> — > 1"0~ is not in J/if> but in if)'. The 
large relative phase between a and a 7 cannot be attributed to a resonance in the s-channel 
of I = 0. Since the 1 = amplitudes receive no net enhancement, a contribution of an 
s-channel resonance, if any, would be a tiny fraction of the whole amplitude pUT3||. Then 
the 1 = amplitudes could not have a large phase close to 90°. It should be pointed out that 
there exist other attempts to explain the relative suppression of if)' — > 1~0~ with different 



dynamical assumptions or intricate dynamical coincidences ||15l . |lq . |i7||18j . The possibility of 
a destructive interference in if)' [jnj, though it is fortuitous, cannot be ruled out in view of 
our finding of the large long-distance final-state interaction in J /if). Whatever the cause of 
the if)' — > 1~CT suppression may be, the group theoretical parametrization of the amplitudes 
remains the same for J /if) — » 1~0~. Though we are unable to choose the solution to this 
if)' — > pir puzzle among the existing models at present, we are confident that this puzzle does 
not interfere with our analysis in this paper. 



IV. IMPLICATION OF THE LARGE RELATIVE PHASE 

In our analysis we have found the first evidence for a large final-state interaction phase 
in a heavy particle decay which is quite different in nature from the common subchannel 
resonant phases. What generates the large relative phase 5 1 between a 7 and a ? It is 
obvious that it must arise from long-distance strong interactions. The short-distance final- 
state interaction phase difference can be evaluated in the quark-gluon picture of perturbative 
QCD. It is of 0(a s /7r) where a s /ir is 0.1 or less. The large phase difference close to 90° found 
in our analysis cannot be produced with the perturbative QCD interaction. The source of 5 7 
must be in the long-distance part of strong interactions, namely, rescattering among hadrons 
in their inelastic energy region. 
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When many channels are open for strong interaction rescattering, the phase of a decay 
amplitude into a physically observed state is determined by the phase shifts of eigenchannels 
of S-matrix and the coupling to them. In the case of the J/ip decay, the decay amplitude 
D(J/ip — > h) into a final state h {e.g., p + n~) is written in the form 

D(J/i> -+h) = a)e iSa Oah, (16) 

a 

where a refers to the eigenchannels of the partial- wave S-matrix of J PC = 1 at the energy 
of the J/ip mass with <5q,'s being their eigenphase shifts. Time-reversal invariance requires 
that X(J/ip — > a)* = X(J/ip — > a) and that O a h be an orthogonal matrix relating h to 
a by (h\ = J2( a \O a h- The partial- wave phase shifts contain much of long-distance physics 
no matter how high the energy is. One indication of substantial long-distance physics in 
the high-energy phase shifts was pointed out by making a partial-wave projection of the 
diffractive scattering amplitude Q. Long-distance physics enters the eigenchannel matrix 
O a h as well. When we compare the three-gluon and one-photon decay amplitudes of J / ip with 
Eq. (|TED , we see no simple relation between their phases in general: The eigenphase factors 
e Sa are summed with the weights X(J/ip — > a) different for J /if) — > ggg and J/ip — > 7 — > qq, 
leading to two phases practically unrelated to each other. In this picture the final-state 
interaction phases of the J/ip decay are generally not determined by short-distance physics 
alone even though pertubative QCD applies to the inclusive J/tj) decays. The analysis of 
this paper indicates that long-distance physics can be far more important in the exclusive 
decays. 

The conclusion of our amplitude analysis, if it is sustained, has a significant implication 
in a wide range of phenomena. For instance, when we evaluate the baryon asymmetry in 
the early Universe from CP-violating particle decays, we compute only the short- distance 
contribution of final-state interactions. Such a calculation makes sense only as an order-of- 
magnitude estimate at best. In the case of the baryon asymmetry we may not ask for a high 
precision after all. However, in the B-meson decay where knowledge of much higher precision 
will be needed for final-state interaction phases, we shall have to know the long-distance final- 
state interaction phases above the inelastic thresholds. It is nearly an impossible task to 
either compute them theoretically or extract them from scattering data. If this is the case, 
the parameters of the fundamental interactions can be extracted only from those data which 
are free from complications due to the final-state interaction. It will not be an easy task to 
look for meaningful physics in the rest of data. 

To conclude this paper, we should emphasize that numerical conclusion of our analysis 
replies on the current data listed in RPP, not only their central values but also the exper- 
imental uncertainties. We cannot rule out the possibility that a future change in the data 
may upset our conclusion, i.e., the need of a large rescattering phase. For this reason, high 
precision measurement of the J/ip decay branchings, particularly for pir, K*K and ujtt , will 
be very important to our understanding of the final-state interactions in general. 
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TABLES 



TABLE I. The SU(3) parametrization of the 1~CT decay amplitudes. Listed are the coefficients 
of the amplitudes written on the top of each column. For instance, the coefficient of a 7 for K*°K 
is —2 x (1/3). The 77 — 77' mixing and the p — uj mixing are introduced as explained in the text 



Decay modes 


a 


e 


3°7 


1, 

3 7 1 


1, 

3 C 72 


a' 


e' 


f 


p + ir~(= p~Tt + ) 


1 





1 

















A° 


1 





1 

















K*+K-(= K*-K+) 


1 


1 


1 


2 














n m 


1 


1 


-2 


-1 














coir 








3 

















P°V8 








V3 





-Vs 


















v/T/3 


















v/2/3 


v/8/3 


-v/8/3 


-v/8/3 


-v/2/3 











p°m 























v/3/2 


UJT)l 

















v/2/3 





VT76 


4>vi 

















-v/173 


-v/173" 


v/173 


07T° 



























TABLE II. The observed branching fractions and our fits with and without 77' included. No 
phase means a fit assuming that all amplitudes be real. The bottom row lists x 2 f° r eacri fit- All 
numbers are in percent except for x 2 - 



Decay modes 


Observed 
(in percent) 


^9^8 


^8,1 


Best fit 


No phase I 


No phase II 


Best fit 


No phase 


p + ir~(= p~Tt + ) 


0.43± 0.03 


0.44 


0.44 


0.47 


0.43 


0.42 


A 


0.42± 0.05 


0.44 


0.44 


0.47 


0.43 


0.42 


K*+K-{= K*-K + ) 


0.25± 0.02 


0.25 


0.30 


0.25 


0.25 


0.31 




0.21± 0.02 


0.21 


0.14 


0.16 


0.21 


0.16 


con 


0.042±0.006 


0.047 


0.034 


0.33 


0.052 


0.030 


<M° 


< 6.8 x 10~ 4 

















p°?7 


0.0193±0.0023 


0.0174 


0.0104 


0.0186 


0.0140 


0.0144 


u 77 


0.158±0.016 


0.131 


0.129 


0.146 


0.146 


0.150 


4>v 


0.065±0.007 


0.065 


0.062 


0.076 


0.064 


0.058 


P°V' 


0.0105±0.0018 








0.0098 


0.0112 


uj-q' 


0.0168±0.0025 








0.167 


0.0169 


H 


0.033±0.004 








0.033 


0.032 


x 2 




4.8 


57 


21 


8.2 


43 
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